Propagation dynamics of off-axis noncanonical vortices in a collimated Gaussian beam
1. IntroductionThere has been considerable interest on the phase singularities nested in wave fronts, i.e., the optical vortices, which are surrounded by wave fields with helicoidal phase structures and may carry the orbital angular momentum.[1] The study on the optical vortices is important both fundamentally and practically, whilst the potential applications include laser trapping,[2] optical data storage,[3] optical information processing,[4] and so on. Optical vortices are also intriguing because they behave in some degree as “particles” in the host beam, and exhibit rich propagation dynamics, while the paths of the vortices can be twisted, looped, and even knotted.[5] When the vortex trajectory is parallel to the wavevector k of the host beam, the resulted beam has helicoidal wave fronts; the situation would be much complicate when the vortex lines take other forms. Consequently, the dynamics of the vortices have been extensively investigated, especially in nonlinear systems,[6–9] including Kerr-type media,[10] cubic–quintic saturating nonlinearity,[11] and nonlocal nonlinear media.[12,13] Even in the linear regime, the free-space propagation of optical vortices may present unusual dynamical properties. For example, identical vortices are found to orbit each other and diffract in a linear medium.[14] Additional vortices can be created when the off-axis vortex is interacted with host beams with dislocation surfaces, e.g., Laguerre–Gaussian (LG) beams.[15] Passing through a cylindrical lens, the noncanonical evolution of an initially canonical vortex leads to a competition between the angular momentum eigenmodes, and results in the inversion of the topological charge.[16] Gouy phase is applied to explain the rotation of multi off-axis vortices that are generated in the collinear superposition of LG beams.[17] The propagation dynamics of both axial[18] and off-axis[19] noncanonical vortices have been investigated. It is pointed out that an off-axis noncanonical vortex moves across a background Gaussian beam without rotation.[14,19]
This paper focuses on the propagation dynamics of a single off-axis vortex in Gaussian beam with or without cylindrical symmetry. It is found that the vortex swirls a while before it approaches infinite, which is mainly due to the divergence of the host beam. By neglecting the divergence, an elliptical orbit around the beam axis can be derived analytically for the symmetrical beam. On the contrary, the vortex trajectory is much complicate in the background of asymmetrical Gaussian beam due to the inversion of the topological charge. In both cases, the vortices are confined in a rectangular area, which is specified by the noncanonical parameter à and the ellipticity of the host beam. To treat the swirling vortex as a “particle”, we define a particle-like angular momentum Lp of the vortex center relative to the beam axis. It turns out that the Lp is conserved only in the symmetrical host beam. Numerical simulation suggests that in the case of asymmetrical beams, the moment of inertia of the off-axis vortex can be roughly approximated as |Re(Ã) |r2. This work may be useful in the exact control over the off-axis vortex position and find application in areas such as optical manipulation, beam shaping, etc.
2. Vortex propagation in symmetric host beamLet us assume a complex scale optical field E(x,y,z) with a noncanonical vortex embedded in it, which can be cast into
where
F(
x,
y,
z) is the host beam, and
describes a vortex of topological charge one, whose sign is the same as that of the real part of the noncanonical parameter
à =
Ar +i
Ai. When
à = ± 1, this vortex reduces to a canonical, i.e., symmetrical vortex, since
x±i
y =
rexp(±i
θ). In Fig.
1, we show a few examples, where arbitrary units are applied.
Two important expressions can be derived from the paraxial wave equation in normalized units
, which determine the propagation dynamics and the variation of the noncanonical parameter[19]
where
L1 =
∂ /
∂x+i
̶/
∂y, and
L2 =
i̶/
∂x–
∂/
∂y. Consider a symmetrical Gaussian beam propagating in a linear homogeneous medium
where
ρ2 =
x2 +
y2,
w0 is the waist radius, and
Z0 is the Rayleigh distance. In view of Eqs. (
4) and (
5), it shows that
à is constant. Substituting Eq. (
5) into Eq. (
3), we obtain the following equation:
which describes the propagation trajectory of the vortex position (
x1,
y1). For the rest of the paper, the initial position of the vortex is denote as (
x10,
y10). Equation (
5) confirms that if the initial position of the vortex is at the center of the beam, then this vortex would not move during propagation. Figure
2 shows the numerically simulated trajectories of several off-axis vortices, where the parameters are
w0 = 1,
Z0 = 10, and
à = 2. It is clear that these vortices will eventually move to infinity, but the obtained results are slightly different from those in Ref. [
19]. The reference predicted that the vortex leaves the host beam in a straight line, but figure
2 shows that the vortices move in a spiral way. The difference is due to the specific approximation used in Ref. [
19]. Careful investigation shows that the initial direction of motion of the vortices in Fig.
2 is perpendicular to the direction of the sloping field.
[14] If the real part of
à is negative, the projections of the vortex trajectories onto the transverse plane will be a large clockwise rotation. This paper will focus on this spiraling property of the off-axis vortex in the linear homogeneous medium, so the divergence of the host beam will be ignored.
Assume that Z0 ≫ z, so the term 1 + iz/Z0 in Eq. (5) can be omitted, which can be recast into the well-known form
where
wx =
wy =
w0 for a symmetrical Gaussian beam, and
wx ≠
wy for the asymmetrical case. First, we consider the simple symmetrical case, where equation (
4) yields d
Ã/d
z = 0, indicating that the noncanonical parameter is invariant. Then equation (
6) can be split into
where the parameter
![](cpb_28_3_034205/cpb_28_3_034205_ieqn2.gif)
. Consider the simple case of
Ai = 0, it immediately follows that the vortex trajectory in the transverse plane is an ellipse, while the semi-major and semi-minor axes coincide with the
x and
y axes, respectively. Consequently, we can write the elliptical trajectory as
where
a and
b are the semi-major and semi-minor axes, respectively. Combining Eqs. (
8) and (
9), it is not difficult to derive
As can be seen from Eq. (
10), the rotation speed of the vortex is determined by the waist radius of the beam, and the rotation direction depends on the sign of
Ar. And the absolute value of the noncanonical parameter
à is equal to the ratio of the major and the minor axes of the ellipse.
For the more general case of Ai ≠ 0, it can be proved that the vortex trajectory is still an ellipse, whose major axis does not need to overlap with the x or y axis. Set the angle between the a axis and the x axis to be ϕ, then we can use a different set of coordinate system (x′,y′), which is related to the original (x,y) by a counter-clockwise rotation ϕ. In the new coordinates, the trajectory is
where
By combining the two equations above, it is not difficult to derive the vortex trajectory in the original coordinate system as
Thus, the trajectory can be solved analytically as follows:
Now we have analytically solved the propagation trajectory of an off-axis vortex, with arbitrary complex noncanonical parameters. Figure
3 gives a concrete example of the dynamics of the rotating vortex around the beam center, whilst the field distribution of the Gaussian beam varies accordingly at different transverse planes. If the host Gaussian beam is well collimated, i.e.,
Z0 ≫ 0, it is possible to observe such rotation before the vortex leaves the beam.
It is interesting to note that there is a common rule related to the vortex trajectory, which still holds even when the host beam becomes asymmetrical. As can be seen from Fig. 4, the elliptical orbits of the vortices locate in rectangular regions of |x1|< xmax and |y1|< ymax. Now we will calculate the ratio between xmax and ymax via Eq. (13)
Substituting Eq. (
14) into Eq. (
15), and through purely algebraic process, yields a surprisingly simple result
Figure
4 plots the numerical simulated projections of the vortex trajectory onto the transverse plane, where the initial positions of the vortices are set at (0,0.8). By varying the phase
ϕA, and keeping the |
Ã| constant, the elliptical trajectories vary differently, but the values of
xmax/
ymax are invariant.
Because the optical vortex behaves in some degree as a particle, we can define an angular momentum of the rotating vortex core relative to the axis of the host beam. Note that this angular momentum is different from the orbit angular momentum carried by the wavefront of the vortex or the whole beam. So we can define the angular momentum of the particle-like vortex core as
where
ez is the unit vector along the
z-axis. Inserting Eq. (
13) into Eq. (
17) yields
And in view of Eq. (
14), it is clear that the angular momentum is completely determined by
à and the initial position (
x10,
y10). Equation (
18) confirms that
Lvc is invariant for the symmetrical host beam, since
à is constant.
3. Vortex propagation in asymmetric host beamFrom now on, we will consider the case of asymmetrical beam, which is given by Eq. (7) with wx ≠ wy. Equation (4) turns into
which follows
where
![](cpb_28_3_034205/cpb_28_3_034205_ieqn3.gif)
. In the asymmetric host beam,
à rotates in the complex plane according to Eq. (
20), whilst the topological charge has the same sign as the real part of
Ã, so the topological charge switches between positive and negative as the beam propagates, this is known as the effect of topological charge inversion. Assuming
Ã(0) = |
Ã|exp(i
ϕA),
![](cpb_28_3_034205/cpb_28_3_034205_ieqn4.gif)
,
![](cpb_28_3_034205/cpb_28_3_034205_ieqn5.gif)
, we can derive the vortex trajectory in the asymmetrical beam as
where
It is interesting to note that the determinant of the matrix is
Equations (
20)–(
22) describe how the vortex moves as the beam propagates and how the topological charge varies. For illustration, figure
5 simulates the dynamics of an off-axis vortex embedded in an asymmetrical beam, the parameters are (
x10,
y10) = (1,0.5),
wx = 1,
wy = 0.8, and
Ã(0) = 1+i. It is clear that the vortex simultaneously rotates at a constant speed and moves along a complicate trajectory.
There is no need to go to the details of the trajectory, which is too complicate. Instead, let us exam if the vortex trajectory is confined in a finite space in the transverse plane. We can calculate the acceleration of the vortex core, which can be derived as
where
It can be proved that
![](cpb_28_3_034205/cpb_28_3_034205_ieqn6.gif)
, and the case of symmetrical host beam falls into the special case of
kA = 0. Let the same vortex in Fig.
5 propagate for a longer distance 0 <
z < 30. Figure
6 plots the projection of its trajectory on the transverse plane. It is obvious that the trajectory is still confined in a rectangle, which proves that the off-axis vortex would not leave the beam if the divergence of the beam is ignored. Imagine that the vortex is at the corner of the rectangle, then both derivatives d
x1 / d
z and d
y1/d
z should approach zero. Inserting this assumption into Eq. (
21), we can derive
which is verified by the numerical results in Fig.
6. And it is not difficult to see that equation (
16) is only a special case of Eq. (
26) when
wx =
wy.
Finally, we consider the angular momentum of the particle-like vortex center relative to the center of the host beam based on Eq. (17), which is shown in Fig. 7. For comparison, we plot the real part of the noncanonical parameter Re(Ã), whose sign defines the topological charge of this vortex. It is found that Re(Ã) has the same period as Lvc, and the latter is no longer invariant for the asymmetrical host beam.
In order to study the physics behind, we further calculate the product of Re(Ã)× |r|2, it turns out that the curves of Lvc and Re(Ã)×|r|2 should be intimately related. It is interesting since these two curves are obtained via quite different equations. The fact that Lvc in an asymmetric host beam can be both positive and negative shows the rotation direction of the off-axis vortex is variable. Equation (21) connects the velocity v = (vx,vy) with the vortex core position r = (x1,y1). It is not difficult to verify that if vx y1 –vy x1 >0, the vortex rotates clockwise, and vice versa. Inserting Eq. (21) into the criteria, it turns out that the rotation direction is determined again by the sign of Re(Ã), which fits well with the simulated results in Fig. 7. As Re(Ã) is negative, the topological charge is −1, the vortex rotates clockwise. It should be note that this rule holds for both symmetric and asymmetric cases (see Fig. 3). Two hypotheses are available here.
(I) If the vortex is treated as a particle, then the term |Re(Ã)|r2 can be used as a rough approximation of the moment of inertia I = mr2, where the effective mass is approximated as |Re(Ã)|. In a symmetric host beam, the effective mass is constant.
(II) The topological charge inversion is the physics behind the rather complicate propagation property of a noncanonical vortex in an asymmetric host beam, which leads to the frequent variation of rotation direction.
4. ConclusionWe analytically investigate the propagation dynamics of a single off-axis vortex nested in both symmetrical and asymmetrical host beams in details, where the vortex swirls around the beam axis if the divergence is ignored. For the symmetrical case, a simple elliptical orbit in the transverse plane can be defined, whilst the topological charge and the angular momentum Lvc of the rotating vortex center are invariant as the beam propagates. For the asymmetrical background beam, the trajectory of the vortex is rather complicate, since the noncanonical parameter varies simultaneously. The topological charge and the angular momentum Lvc relative to the beam axis are not conserved. But in both cases, the vortex is always confined in a rectangular area and the rotation direction is determined by the sign of the topological charge.